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The investigation of interactions between magnetic nanotubes is complex and often involves sub- 
stantial simplifications. In this letter an analytical expression for the magnetostatic interaction, 
taking into account the geometry of the tubes, has been obtained. This expression allows for the 
definition of a critical vertical separation for relative magnetization between nanotubes and can 
be used for tailoring barcode-type nanostructures with prospective applications such as biological 
separation and transport. 



Since the discovery of carbon nanotubes by Iijima in 
1991, 1 intense attention has been paid to tubular nanos- 
tructures. Because of their geometry, nanotubes offer 
prospective applications in catalysisj^GD sensors,^ and bio- 
logical separation and transport pEHZJ i n particular, mag- 
netic nanotubes are the focus of growing interest due to 
the existence of techniques that lead to the production of 
highly ordered arrays PEll These particles offer an addi- 
tional degree of freedom as compared to nanowires; not 
only can the length, L, and external radius, R, be var- 
ied, but also the internal radius, a. In this way, nan- 
otubes may be suitable for applications in biotechnology, 
where magnetic nanostructures with low density, which 
can float in solutions, become much more useful for in 
vivo applications.^ These tiny magnetic tubes could pro- 
vide an unconventional solution for several research prob- 
lems, and a usef ul vehicle for imaging and drug delivery 
applications 

In such systems changes in thickness are expected to 
strongly affect the mechanism of magnetization reversal 
and, thereby, the overall magnetic behavior. Also inter- 
actions play a fundamental role, modifying the magnetic 
behavior of the particles. Clearly, for the development 
of magnetic devices based on those arrays, knowledge 
of the magnetostatic interaction between the tubes is of 
fundamental importance. But as usual the effects of in- 
terparticle interactions are complicated by the fact that 
the dipolar field felt by each element depends upon the 
magnetization state of all the elements in the array. In a 
previous work by Lee et. a/, pi multisegmented metallic 
nanotubes with a bimetallic stacking configuration along 
the nanotube axes were prepared and investigated. These 
particles exhibit different magnetic behaviors, which en- 
courage a study about the magnetostatic interactions be- 
tween the stacking. Due to the very narrow hysteresis 
loops that are obtained, the influence of the interactions 
is not easily identifiable from magnetization curves, and 
then a theoretical study can shine light on this problem. 

The purpose of this letter is to develop an analytical 
model for the full long-range magnetostatic interaction 
between two nanotubes exploring the possibility of vary- 
ing the magnetic coupling as a function of the tubes posi- 
tion. The geometry of the tubes is characterized by their 
external and internal radii, R and a, respectively, and 



length L. It is convenient to define the ratio /3 = a/R, so 
that (3 = represents a solid cylinder (wire) and (3 close 
to 1 corresponds to a tube with very thin walls. The 
separation between the tubes is written in terms of the 
inter-axial distance, and the vertical separation, s, as 
depicted in Fig. 1. Our model goes beyond the dipole- 
dipole approximation and lead us to obtain an analytical 
expression for the interaction in which the lengths and 
radii of the tubes are taken into account. We focus on 
the stability of parallel and antiparallel magnetization 
alignment in pairs of interacting tubes, as a function of 
the distance between them, in order to gain insight on the 
understanding of the role of interactions on barcode-type 
nanotubes.™^ 
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FIG. 1. (a) Geometric parameters used for the individual 
tube description, (b) Relative position of interacting tubes: 
d is the inter-axial distance and s the vertical separation. 

We adopt a simplified description of the system, in 
which the discrete distribution of magnetic moments is 
replaced with a continuous one characterized by a slowly 
varying magnetization M(r). The total magnetization 
can be written as M (r) = Mi (r) + M 2 (r), where Mi (r) 
and M2 (r) are the magnetization of tubes 1 and 2, re- 
spectively. In this case, the magnetostatic potential U (r) 
splits up into two components, U\ (r) and U2 (r), asso- 
ciated with the magnetization of each individual tube. 
Then, the magnetostatic energy of two interacting mag- 
netic tubes may be written in terms of their magnetiza- 
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tions and the fields generated by each one. The general 
expression, after using the reciprocity theorem, is E d = 
EL + El K + E int . The El u =<fjM i (r). W; (r) dV, 
terms correspond to the self-energy of the i-th tube, and 
E^t = f M 2 (r) • VE/i (r) dV 2 is the interaction energy 
between tubes, which is the focus of this letter. 

In order to proceed, we first need to calculate the mag- 
netostatic potential U (r) of a single tube. However, it is 
necessary to specify the functional form of the magnetiza- 
tion for each nanotube. Due to their geometry and in or- 
der to reduce the stray field, tubes present what is called 
a flower configuration.^ In this case, the magnetization 
is mostly homogeneous, and spreads outwards near the 
ends. This non-homogeneity produces a decrease of the 
interaction energy felt by one tube due to the other. For 
L^> R this decrease is small and can be neglected to sim- 
plify the calculations. Thus, we consider tubes with an 
axial magnetization defined by (r) = Moa^z, where 
Mo is the saturation magnetization of each nanotube, z 
is the unit vector parallel to the axis of the nanotube 
and Gi takes the values ±1, allowing the magnetization 
of tube i to point up (c^ = +1) or down (c^ = —1) along 
z. The magnetostatic potential produced by the tube 1 
with magnetization Mi(r) is given by volume and surface 
contributions and can be written as 
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Note that the first term on the right-hand side of Eq. 
vanishes because the magnetization field is constant. Fur- 
thermore, the surface integral in Eq. ([!]) has contributions 
only from the upper and lower ends of the tube, located 
at z = L/2 and z = —L/2, respectively. Due to the sym- 
metry of the problem, we have calculated the integral 
in ([l]) using a kernel in cylindrical coordinates.^ After 
some manipulation, the integral expression for the scalar 
potential can be written as 
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where J n (x) is the n-th order Bessel function of first kind. 
Figure 2 shows the surface plot of U\ for a thin- walled 
tube (/? = 0.8) and a wire (/? = 0). From this figure it 
is clear that one should expect distinct behaviors from 
these two nanoelements in the region near the ends. 

Now it is possible to calculate the magnetostatic inter- 
action energy between two identical nanotubes using the 
magnetostatic field experienced by one of the tubes due 
to the other. The final result reads 
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FIG. 2. Magnetostatic potential U (r) according to equation 
(2) produced by a tube with L = 500 nm, R = 50 nm, and 
(a) f3 = 0.8 (thin- walled tube) and (b) = (wire). The 
color scale is chosen such that higher absolute values of U are 
represented by darker shades. 



Here E [nt is the interaction energy in units of [IqMqV 
i.e. E int = E int /n MgV, where V = ttR 2 L (l - j3 2 ) is the 
volume of a tube. Equation (3) has been previ ously ob- 
tained for nanowires (/? = 0) for s = nm! 18 * 19 * The gen- 
eral expression for the interaction energy between tubes 
with axial magnetization, given by equation (3), can only 
be solved numerically. However, tubes that motivated 
this workPESl satisfy Rj L = a <C 1, in which case one 
can use that J\ (ax) « ax/2. With this approximation, 
Eq. (3) can be written in a very simple form as 



E int = 
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The simplicity of Eq. Q makes it an excellent tool for 
the understanding of interactions between those nanoele- 
ments. Figure 3 illustrates the interaction energy, ob- 
tained from Eq. (4), between two identical nanotubes 
with parallel axial magnetization as a function of 2R/d. 
When the two tubes are in contact, 2R/d = 1; when 
they are infinitely separated, 2R/d = 0. Differences be- 
tween this expansion and the full expression in Eq. (3) 
are less than 3% for any L/R. As an illustration, when 
we consider two Ni nanotubes (/? = 0.8) with L = 500 
nm, i? = 50 nm, s = 100 nm and d = 100 nm, we obtain 
Eint = 13.15 eV from Eq. (3), and Emt = 13.75 eV from 
Eq. (4). 

It is interesting to analyze the behavior of the in- 
teraction energy given by Eq. Q as the inter-axis dis- 
tance, is kept fixed and the vertical separation, 8, is 
varied. We define the critical separation so such that 
E int (L,d,so) = 0. For tubes with L = 500 nm and 
d = 100 nm, so = 317 nm, independently of other pa- 
rameters, as depicted in Fig. (4). Since lateral positions 
(small s) favor antiparallel magnetization alignment, the 
interaction energy is positive for s < sq = 317 nm. For 
s > so = 317 nm the interaction between tubes with par- 
allel magnetization is atractive, the strongest attraction 
appearing for s ~ 500 nm. 



3 



0.02 



0.00 



F ■ 


1 1 1 ■ X 




— p-ooo 






p-0.40 






p-0.80 


/ s 




M99 


/ s 
/ s 






/ s _ 




/ s 
/ s 




/ s 




/ s 




' s . • ' 











0.0 0.2 0.4 0.6 0.8 1.0 

2R/d 



FIG. 3. (Color online) Interaction energy between two identi- 
cal nanotubes, as given by Eq. (4). The tubes have L = 500 
nm, R = 50 nm, different values of /3 and parallel magneti- 
zation defined by <j\ — 02- The vertical position kept fixed, 
s = nm, and the inter-axis distance d is varied as a function 
of 2R/d. 



be seen in Fig. (5) For L > 2.5 /im, which corresponds 
to values usually found in nanotubes, we observe a linear 
dependence of the form sq = 0.62L, almost independent 
of d. 

In conclusion, by expanding the general expression for 
the magnetostatic interaction energy between tubes, and 
keeping the first-order term we have obtained an expres- 
sion that can be easily used to calculate the magneto- 
static interaction between tubes. In particular, we have 
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FIG. 4. (Color online) Interaction energy of two identical 
nanotubes, as given by Eq. (4). The tubes have L = 500 nm, 
R = 50 nm, different values of ft and parallel magnetization 
defined by <j\ =02- The inter-axis distance kept fixed, d — 
100 nm, and the vertical separation s is varied. 



The dependence of sq on L for different values of d can 



FIG. 5. (Color online) Dependence of the critical vertical 
separation so on d and L for two identical tubes with arbitrary 
values of f3 and R. For s < so antiparallel alignment of the 
magnetization is favored. 



investigated the relative position of the tubes for which 
the ferromagnetic and antiferromagnetic configurations 
are of lowest energy. For the tubes usually found in the 
literature, L > 2.5 /im, we observe a linear dependence 
of the critical vertical distance of the form so = 0.62L. 
Our results are intended to provide guidelines for the pro- 
duction of barcode-type nanostructures with prospective 
applications such as biological separation and transport. 
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